Exchange and the Coulomb blockade: Peak height statistics in quantum dots 
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We study the effect of the exchange interaction on the Coulomb blockade peak height statistics in chaotic 
quantum dots. Because exchange reduces the level repulsion in the many body spectrum, it strongly affects 
the fluctuations of the peak conductance at finite temperature. We find that including exchange substantially 
improves the description of the experimental data. Moreover, it provides further evidence of the presence of 
high spin states (S > 1) in such systems. 

PACS numbers: 73.23.Hk, 73.40.Gk, 73.63.Kv 
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Transport properties of quantum dots (QDs) are strongly 
affected by interactions. A paradigmatic example is the 
Coulomb blockade (CB) of electron tunnelingpi^ It occurs 
at low temperature when the thermal energy is smaller than 
the charging energy {Ec) required to add an electron to the 
QD. The conductance is then blockaded, being restored only 
at specific values of the gate voltage V g of a nearby gate ca- 
pacitively coupled to the QD. This leads to a series of sharp 
CB conductance peaks. The statistical properties of these 
peaks reflect the mesoscopic fluctuations of the (many-body) 
spectrum and wave-functions of the QD. In particular, they 
encode information about the non-trivial spin statistics that 
arises from the presence of exchange, the most important part 
of the residual interaction iM 

Both the peak height distribution (PHD) and the peak spac- 
ing distribution have been studied in detail (see Ref. |2| for a 
review), though the latter has received most of the attention. 
This is mainly due to the fact that early measurements^ of 
the PHD were found to be in good agreement with the con- 
stant interaction (CI) model, 7 8 contrary to the case of the peak 
spacing. In this simple model of CB, the e-e interaction is as- 
sumed constant (given by Ec) and the fluctuations of the sin- 
gle particle properties described by random matrix theory — 
appropriate for chaotic (or diffusive) QDs. 

A later experiment 9 showed, however, significant devia- 
tions from the CI model. Namely, the peak height fluctuations 
were found to be smaller than expected in the entire experi- 
mental temperature range (k B T = 0. 1 — 2 A, with A the single 
particle mean level spacing). The authors attributed the dis- 
crepancy at high T to the presence of dephasing but the ori- 
gin of the low T behavior was unclear — a recent calculation 10 
showed that dephasing by itself cannot account for the high T 
data either. It was then proposed 11 that the observed reduc- 
tion of the fluctuations at low T is related to the presence of 
spin-orbit coupling. Though the result is in good agreement 
with the low T data, the assumed strength of the spin-orbit in- 
teraction requires some test — for instance, experimental data 
in Ref. 12 suggest that in small QDs (as used in Ref. 0) the 
effect of spin-orbit is rather weak. 

So far, however, exchange — the main interaction effect — 
has not been considered. This contrasts with the case of 
the peak spacing distribution where it was shown to be cru- 
cial, as it leads to the appearance of non-trivial spin states 
in the QD (S > l) 13 i 14 i 15 i 16 and enhances the effect of finite 
temperature»i2ii£ 



In this work we show that exchange also strongly affects the 
PHD. Furthermore, it accounts for most of the disagreement 
with the low T data of Ref. 9, without including dephasing or 
spin-orbit coupling. Such agreement is yet another indication 
of the presence of high spin states in QDs. 

A chaotic QD containing a large number of electrons can 
be described by the following HamiltonianiiA^ 



-E c (n-Af) 2 -J s S 2 



(1) 



where {e Q }, the single particle energies, are described by ran- 
dom matrix theory, J\f=C g V g /e describes the capacitive cou- 
pling to the control gate, C g is the dot-gate capacitance, S 
is the total spin operator, and J$ is the exchange constant. 
The second and third terms in Eq. Q describe the QD charg- 
ing energy and the exchange interaction, respectively. We as- 
sume elastic transport and broken time-reversal symmetry — 
therefore we use the Gaussian Unitary Ensemble (GUE) — 
unless otherwise stated. 

In the regime T < k B T,A < E c , where T is the total 
width of a level in the QD, the conductance near the CB peak 
corresponding to the N— 1 — > N transition is given byi&i2i22i2I 
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where x = 2Ec[N— (N— i)] and x* % is defined below. Here, 



(i) T„ is the partial width of the single-particle level a due 
to tunneling to the left (right) lead, (ii) {| v I'j v )} are the eigen- 
states of the QD with N electrons, (iii) F eq (j) is the canonical 
probability that the eigenstate j is occupied, and (iv) 



f(yii-x)f(yji+x) 



(3) 



- exp[x/k B T}) 1 and = 
Notice that gji(x) reaches its 
0. Let us denote by {Ej} the 

eigenvalues of H without the charging energy term. Then, the 
contribution of the transition i— *j to the conductance reaches 



with f(x) = (1 
k B Tln[F eq (j)/F eq (i)}. 
maximum value at x = 
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FIG. 1 : Coulomb blockade peak height distribution at finite temper- 
ature. The theoretical distributions (lines) were obtained by numer- 
ically solving Eq. J2j for different values of the exchange constant 
Ms). The histograms correspond to the experimental data in Ref. 
9. Notice that good agreement is obtained only after the addition of 
exchange. 



its maximum when x=x* H = Ef—E^ 1 +y, l /2-E F . Note the 
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shift in the peak position due to finite temperature 

Since both F eq and the overlap I |c+ J^f -1 ) I depend 
on the spin of the two states involved in the transition, it is 
clear that this maximum is spin-dependent>I2iSSiSi21 

The eigenstates of the Hamiltonian Q can be written as 
l*f ) = IK>- S', S' z ,k'), with n' a = n' aA +n' ad the occupa- 
tion of the a-th single particle state. 42 Since there is no en- 
ergy dependence on S' z — Zeeman energy is neglected — the 
sum over S z , S' z and a in Eq. (EJ can be easily carried out. 
One then gets 1 ^ 



E 



N 

St,k< 



N-l\ 
hi 



'(2S' + l)N k ifn Q = 
(2S+l)N k < ifn a = l 



(4) 

where n a refers to the state with N— 1 particles, Nk (Nk>) is 
the degeneracy associated with fc (fc') 43 an d | *S" — /S" | = \. 

In the special case of very low T, when only the ground 
state is relevant, we have F eq (j) ~ 1/(25" + !) , F eq (i) — 



1/(25+1), and N k ,=N k = l, 
given by 



so that the peak conductance is 
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where still |5"— S\ — \- Because the probability for the S^S' 
transition depends on both the interaction strength ( Js) and 
the statistics of the single particle spectrum, it turns out that 
the PHD depends, even at low T, on the fluctuations of both 
the wave-functions (Ps) and the many-body spectrum. 

At finite temperature many transitions contribute to the con- 
ductance. We calculate the PHD by finding the maximum of 
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FIG. 2: Ratio of the root mean square to the mean of the conductance 
peak as a function of temperature. Different curves correspond to 
different values of exchange Js ~0, 0.3A, and 0.4A (dotted, solid, 
and dashed lines, respectively). Symbols correspond to the data in 
Ref. 0. Notice the strong reduction of the fluctuation introduced by 
the exchange interaction at low temperature. The inset shows the 
decay of o-(G P cak)/{G P i:ak) as a function of Js at ksT — 0.35A. 
At high T the experimental data show a stronger suppression of the 
fluctuations than predicted in the strong inelastic regime for Js = 
0.4A (dot-dashed line). 



Eq. numerically. According to random matrix theory, we 
describe the fluctuation of the widths with q = L, R by 
the Porter-Thomas distribution P(r«) = f exp(-r« /f ), 
appropriate for the GUE. We keep states within an energy win- 
dow of max(2A, 6/cb? 1 ) around the ground state. 

Figure[T]shows the PHD for k B T = 0.1 A, 0.5A and J s = 0, 
0.3A, 0.4A. The histograms correspond to the experimental 
data in Ref. 9. At low temperature the agreement is very 
good for the non-zero values of the exchange constant but 
clearly not for Js — (CI model). This is indirect evidence 
for the presence of high spin states in QDs — we should point 
out though that spin-orbit leads to a similar effect if it is strong 
enough. 1 1 At higher temperature the presence of exchange im- 
proves the fit but not enough to fully account for the observed 
distribution. Nevertheless, it is clear that exchange substan- 
tially modifies the PHD and thus cannot be ignored. 

Since an accurate measurement of the full PHD is quite de- 
manding, a detailed comparison with theory is difficult. In- 
stead, it is usually more convenient to look at the first few 
moments. In Ref. 9, the ratio of the root mean square and the 
mean value of the conductance peak height, cr(G pe ak)/ (G pea k), 
was measured as a function of temperature. It was found to be 
smaller than the value predicted by the CI model. Figure [2] 
compares the experimental data with our results for different 
values of Js- The inset shows the Js -dependence at a fixed 
temperature (fc B T = 0.35A). 

The improvement introduced by the exchange interaction 
at low temperature is evident. This can be easily under- 
stood as follows: because exchange reduces the level repul- 
sion in the many body spectrum, the number of levels that 
contribute to the conductance at a given temperature increases 
and therefore the fluctuations are reduced — in fact, the im- 
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FIG. 3: Same as Figure |2| but for zero magnetic field (GOE). The 
GUE result for Js = is included for comparison (dot-dot-dash line). 



portance of the interplay between temperature and exchange 
was pointed out in Ref. [l7l This is similar to the mecha- 
nism discussed in Ref. [IT] where spin-orbit coupling reduces 
the repulsion in the single particle spectrum. At very low 
temperature, Eq. l|5} predicts an enhancement of the fluctu- 
ations compared to the CI model result, er(Gpeak)/(G pe ak) = 
K A S',s)/( A S',s) 2 x 9/5 -l] 1 / 2 > 2/VE. However, this en- 
hancement is negligible since (A|, s ) / (Xs'.s) 2 — 1 f° r typ- 
ical values of the interaction. We also checked that the ef- 
fect on Fig. |3 of non-universal corrections 3 to the Hamilto- 
nian Q is very small, consistent with the results for the peak 
spacing distribution where they mainly affect the shape of the 
distributioni 17 ! 18 ! 26 

Our model shows a significant deviation from the experi- 
mental data for k%T> 0.5 A. It is tempting to attribute this to 
the presence of inelastic processes, which we have not taken 
into account so far. Notice that after exchange is included, de- 
phasing only needs to account for a (small) fraction of the re- 
duction of the fluctuations. Accounting for an arbitrary inelas- 
tic rate requires solving a master equation for the transition 
probabilities, which is not a simple task. 19 ' 27 Instead, we cal- 
culate the peak conductance in the strong inelastic regime, 28 
where electrons inside the QD are assumed to be in thermal 
equilibrium. In practice, this means that instead of calculating 
the thermal average of r^r^/(T^+r^) — see Eq. (0 — , one 
first calculates the thermal average of the couplings {T'} and 
then the ratio, (T L ) (T R ) /(T L +T R ). This is expected to be 
a lower bound for the experimental data. Figure shows that 
this is not the case. The reason for this is unclear. 

A direct comparison between our results and the experi- 
mental data seems to suggest that Js — 0.4A (see inset in 
Figure |3- This value of the interaction is bigger than what 
one would estimate from an RPA calculation of the screened 
potential for the experimental density (Js — 0.3A)M& At 
present it is not clear to us whether corrections beyond RPA 
are required, 29 - 30 or if other effects such as spin-orbit— need 
to be included. 

It is interesting, then, to consider ways to distinguish be- 
tween the exchange and spin-orbit scenarios. In the latter 



FIG. 4: Relative change in the mean conductance upon breaking 
time-reversal symmetry [a, see Eq. J4j] as a function of tempera- 
ture. Different curves correspond to different values of the exchange 
constant: (dotted), 0.3A (solid), and 0.4A (dashed). Notice that a 
is increased by exchange and, in particular, that it is bigger than 0.25 
for T->0 when J s > 0. 



case, it was predicted 11 that an in-plane magnetic field will 
restore the level repulsion and increase the fluctuations. The 
behavior is the opposite in the absence of a perpendicular 
magnetic field, where the system will evolve (as the in-plane 
field is increased) from the orthogonal to the unitary ensem- 
ble (see Figure[3}- 1 1 The same general trend is also valid in the 
exchange-dominated scenario because a large enough parallel 
field will always drive the system towards the strong spin-orbit 
regime. 31 Therefore, the only difference between the two sce- 
narios would be the specific dependence of cr(G pe ak) / (G pe ak) 
on the magnitude of the parallel field. 

A sharper distinction between the two scenarios could be 
made by changing the electron density n e in the QD. On 
the one hand, an increase of n e will slightly decrease Js 
and therefore enhance the fluctuations (see Figures|2]and|3jl. 
On the other hand, increasing n e increases the magnitude of 
the spin-orbit couplin g 12 ! 32 and so decreases the conductance 
fluctuations — according to Ref. [Tj, we assume that the initial 
value of the spin-orbit coupling is smaller than the one re- 
quired to complete the crossover described in Ref. 1 1 . There- 
fore, the sign of 9[er(G P eak)/ (G pe ^}]/dn e at low T is different 
for each scenario — a sharp distinction. 

It is worth pointing out that an independent measure- 
ment of Js — obtained, for instance, by measuring the spin 
distributior. 3 ^^ 3 ^ 6 . — would be a direct way to rule out or 
quantify the effect of spin-orbit coupling and exchange. 

Finally, we study the effect of exchange on the relative 
change of the mean value of the conductance upon breaking 
time-reversal symmetry, 



a=l- 
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Very recently, the temperature dependence of this quantity 
was measured by Folk et al.M who used it to estimate the 
dephasing time (r^) in weakly coupled QDs. 38 This estimate 
is based on the deviation of the experimental value of a from 
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the one expected within the CI model in the elastic regime, 
namely, a ~ 0.25 for both k B T < A and k B T > AX For 
intermediate temperatures a is slightly smaller than 0.25 be- 
cause of the difference in the spectral fluctuations of the two 
ensembles. 10 - 40 

Figure |4] illustrates the temperature dependence of a for 
three different values of Jg. Notice that a increases monoton- 
ically with Jg in the full temperature range. That is, exchange 
enhances a. In particular, the bound a < 0.25 is no longer 
valid at low temperature. 18 This can be verified using Eq. @, 

3 <^5)goe >025 (g) 
4 (As',s)gue 

This enhancement originates in the difference between the 
spin distributions in the two ensembles — high spin states are 
more likely to occur in the GOE. Also note that the effect of 
the spectral fluctuations is substantially reduced. As a result, 
a remains closer to 0.25 than in the CI model. Unfortunately, 
because of large statistical errors in the experimental data, 37 



this difference could not be discerned. 

Since inelastic processes reduce apiS 8 the bigger the dif- 
ference between the elastic result and the measured value of a 
the shorter t$. Then, it is clear that neglecting exchange over- 
estimates T0. Though we do not expect the difference to be 
large, it could be relevant when deciding whether or not there 
is a saturation of t^. 

In summary, we have studied the effect of the exchange in- 
teraction on the peak height distribution. We found that it 
strongly affects the distribution at finite temperature. Our re- 
sults indicate that the experimental data present signatures of 
high spin states and suggest that J$ ^ 0.3A for the dots of 
Ref.g 

We appreciate helpful discussions with K. Held, E. Eisen- 
berg, P. Brouwer, and B. L. Altshuler. GU acknowledges par- 
tial support from CONICET (Argentina). This work was sup- 
ported in part by the NSF (DMR-0103003). After submission 
of this paper, we learned of Ref. 41 which contains some 
overlapping material. 



1 L. P. Kouwenhoven, C. M. Marcus, P. L. McEuen, S. Tarucha, 
R. M. Westervelt, and N. S. Wingreen, in Mesoscopic Elec- 
tron Transport, edited by L. L. Sohn, L. P Kouwenhoven, and 
G. Schon (Kluwer, New York, 1997), pp. 105-214. 

2 Y. Alhassid, Rev. Mod. Phys. 72, 895 (2000). 

3 I. L. Aleiner, P W. Brouwer, and L. I. Glazman, Phys. Rep. 358, 
309 (2002). 

4 I. L. Aleiner and L. I. Glazman, Phys. Rev. B 57, 9608 (1998). 

5 A. M. Chang, H. U. Baranger, L. N. Pfeiffer, K. W. West, and T. Y. 
Chang, Phys. Rev. Lett. 76, 1695 (1996). 

6 J. A. Folk, S. R. Patel, S. F. Godijn, A. G. Huibers, S. M. Cronen- 
wett, C. M. Marcus, K. Campman, and A. C. Gossard, Phys. Rev. 
Lett. 76, 1699 (1996). 

7 R. A. Jalabert, A. D. Stone, and Y. Alhassid, Phys. Rev. Lett. 68, 
3468 (1992). 

8 Y. Alhassid and A. Wobst, Phys. Rev. B 65, 041304 (2002), also 
found good agreement using a model of spinless interacting parti- 
cles. 

9 S. R. Patel, D. R. Stewart, C. M. Marcus, M. Gokcedag, Y. Alhas- 
sid, A. D. Stone, C. I. Duruoz, and J. S. Harris, Phys. Rev. Lett. 
81, 5900 (1998). 

10 T. Rupp, Y. Alhassid, and S. Malhotra, Phys. Rev. B 65, 193304 
(2002). 

11 K. Held, E. Eisenberg, and B. L. Altshuler (2002), cond- 
mat/0208177. 

12 D. M. Ziimbuhl, J. B. Miller, C. M. Marcus, K. Campman, and 
A. C. Gossard, Phys. Rev. Lett. 89, 276803 (2002). 

13 P. W. Brouwer, Y. Oreg, and B. I. Halperin, Phys. Rev. B 60, 
R13977 (1999). 

14 H. U. Baranger, D. Ullmo, and L. I. Glazman, Phys. Rev. B 61, 
R2425 (2000). 

15 I. L. Kurland, I. L. Aleiner, and B. L. Altshuler, Phys. Rev. B 62, 
14886 (2000). 

16 D. Ullmo and H. U. Baranger, Phys. Rev. B 64, 245324 (2001). 

17 G. Usaj and H. U. Baranger, Phys. Rev. B 64, 201319(R) (2001). 

18 G. Usaj and H. U. Baranger, Phys. Rev. B 66, 155333 (2002). 

19 C. W. J. Beenakker, Phys. Rev. B 44, 1646 (1991). 



20 Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 (1992). 

21 H. M. Pastawski, Phys. Rev. B 46, 4053 (1992). 

22 L. I. Glazman and K. A. Matveev, Pis'ma Zh. Eksp. Teor. Fiz. 48, 
403 (1988), [JETP Lett. 48, 445 (1988)]. 

23 H. Akera, Phys. Rev. B 59, 9802 (1999). 

24 D. Weinmann, W. Hausler, and B. Kramer, Phys. Rev. Lett. 74, 
984 (1995). 

25 R. Pauncz, Spin Eigenf unctions: Construction and Use (Plenum 
Press, New York, 1979). 

26 Y. Alhassid and S. Malhotra, Phys. Rev. B 66, 245313 (2002). 

27 E. Eisenberg, K. Held, and B. L. Altshuler, Phys. Rev. Lett. 88, 
136801 (2002). 

28 C. W. J. Beenakker, H. Schomerus, and P. G. Silvestrov, Phys. 
Rev. B 64, 033307 (2001). 

29 H. Jiang, H. U. Baranger, and W. Yang, Phys. Rev. Lett. 90, 
026806 (2003). 

30 Y. Oreg, P. W. Brouwer, X. Waintal, and B. I. Halperin (2001), 
and references therein, cond-mat/0 109541. 

31 J. A. Folk, S. R. Patel, K. M. Birnbaum, C. M. Marcus, C. I. 
Duruoz, and J. S. Harris, Phys. Rev. Lett. 86, 2102 (2001). 

32 J. B. Miller, D. M. Ziimbuhl, C. M. Marcus, Y. B. Lyanda-Geller, 
D. Goldhaber-Gordon, K. Campman, and A. C. Gossard (2002), 
cond-mat/0206375. 

33 P. Jacquod and A. D. Stone, Phys. Rev. Lett. 84, 3938 (2000). 

34 I. L. Kurland, R. Berkovits, and B. L. Altshuler, Phys. Rev. Lett. 
86, 3380 (2001). 

35 J. A. Folk, C. M. Marcus, R. Berkovits, I. L.Kurland, I. LAleiner, 
and B. L. Altshuler, Physica Scripta T90, 26 (2001). 

36 S. Lindemann, T. Ihn, T. Heintzel, W. Zwerger, K. Ensslin, 
K. Maranowski, and A. C. Gossard, Phys. Rev. B 66, 195314 
(2002). 

37 J. A. Folk, C. M. Marcus, and J. S. Harris, Phys. Rev. Lett. 87, 
206802 (2001). 

38 B. L. Altshuler, Y. Gefen, A. Kamenev, and L. S. Levitov, Phys. 
Rev. Lett. 78, 2803 (1997). 

39 Y. Alhassid, Phys. Rev. B 58, 13383 (1998). 

40 K. Held, E. Eisenberg, and B. L. Altshuler, Phys. Rev. B 66, 



5 



033308 (2002). 

Y. Alhassid and T. Rupp (2002), cond-mat/0212126. 

An extra quantum number k' is required when the number of 

singly occupied levels is bigger than 2 (see Ref. l25h . 



If Ns is the number of singly occupied levels, then Nk = 

( N S \ _ ( n s \ 



